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SUMMARY 

The weight function method was used to determine the effective stress 
intensity factor (AK ef f) and the crack opening profile for a fatigue tested 
composite which exhibited fiber bridging. The bridging mechanism was modeled 
using two approaches; the crack closure approach and the shear lag approach. 
The numerically determined AK e ff values from both methods were compared and 
correlated with the experimentally obtained crack growth rates for SiC/Ti-15- 3 
[0] 8 oriented composites. The near crack tip opening profile was also deter- 
mined for both methods and compared with the experimentally obtained 
measurements . 


INTRODUCTION 

The ability to predict fatigue crack growth behavior of metal matrix com- 
posites is of particular interest due to the presence of many cracklike 
defects which tend to be present in these materials. Experimental studies on 
a number of these composites (refs. 1 to 3) have shown that cracks tend to 
propagate in the metallic matrix leaving behind unbroken fibers which bridge 
the two cracked surfaces. 

A number of researchers (refs. 3 to 6) have attempted to model crack 
bridging. The crack bridging model that has received the most attention for 
composite materials is the so-called shear lag model (refs. 3 to 5). This 
model is based on the relative sliding between the fiber and the matrix in the 
region where the interface shear stresses exceed the frictional shear forces. 
This model was developed for brittle matrix fiber composites which have purely 
frictional bonding between the fibers and matrix. This model was also applied 
to metal matrix composite systems (ref. 3). Another type of bridging model 
which may be applicable for these types of composites is based on the applica- 
tion of a closure pressure in the bridged zone proportional to the load car- 
ried by the bridging fibers (henceforth this model will be termed crack closure 
model). This methodology had been applied previously to model the effect of 
unbroken ligaments in a fast fracture test of steels (ref. 6). 



In order to decide whether these types of models can represent the actual 
fatigue behavior, an experimental FCG test program was conducted. Testing was 
performed using a specially designed loading stage mounted inside an SEM. 
this apparatus permits an accurate measurement of the crack opening displace- 
men s, measurement of crack growth rates and identification of active fatigue 
fai lure mechanisms . 6 

Both types of models were used to predict the fatigue crack growth behav- 
ior and the near tip crack displacements. Analytical predictions were compared 
to the test data. r 


EXPERIMENTAL 

Fatigue crack growth tests were conducted on two single edge notch (SEN) 
metal matrix composite specimens. The test material was a SiC (SCS-6) fiber 
reinforced composite with a Ti-l5V-3Cr-3Al-3Sn (Ti-15-3) matrix. The testing 
was performed on unidi rect ional ly oriented specimens with the fiber direction 
being parallel to the loading axis and perpendicular to the starter notch Com- 
panion testing was performed on compact tension (CT) specimens machined out of 

laminated, unreinforced matrix material which was processed in a similar manner 
as the compos 1 te . 


RESULTS 

The fatigue crack growth behavior of the SEN composite specimens consisted 
mostly of stable crack growth perpendicular to the loading axis. Crack growth 
f 3 K con ^ in ® dto t | le matrix ligaments only, leaving behind in its wake unbroken 
fibers bridging the two cracked surfaces (ref. 2). A plot of the crack length 
versus number of cycles (fig. l) revealed a continuous decrease in the fatigue 

C I a ^k gr0W J h rate Wlth increasin § cra ck length. After growth around four row 
ot fibers (approximately 1 mm of growth) the crack growth arrested. 

i ** A j 3 fir f t attempt ’ the fatigue crack growth rate of the composite was 
p otted as a function of the stress intensity factor range (AK ann ) and is 
shown in figure 2 The stress intensity range was calculated by P considering 
only the applied load and the crack length and without taking into considera- 
tion the fiber bridging of the cracked surfaces. The data show crack growth 
ra e ecreasing w i t h an increasing AK a np. Also shown in the figure is the 
fatigue crack growth data of the unreinforced matrix material The crack 
growth rates for the composite are at least a factor of 10 3 lower than those 
ot the unreinforced composite for the same AK ann 

app • 


.ANALYSIS 

Both of the models evaluated use fracture mechanics as the basis for cal- 
culating the effective crack driving force, AK eff . The effective stress inten- 
S1 V actors and the crack opening displacements were determined using the 
weight function method. The weight function used is based on the Bueckner 
(ret. 7) formulation for the stress intensity factor calculation of a single 
edge notch in a finite geometry specimen (fig. 3). The stress intensity fac- 
tor, S1F, tor an unbridged crack of length a is given by: 


2 



dx 


(i) 


K(a) 



p(x)H(a,xJ 
ya - x 


where H(a,x) = (1 + m^a - x)/a + m 2 (a - x) 2 /a 2 ) and where m x and m 2 are 
given by: 

mi = 0.6147 + 17.1844 (a/w) 2 + 8.7822 (a/w) 6 
m 2 = 0.2502 + 3.2889 (a/w) 2 +■ 70.0444 (a/w) 6 
The unbridged displacement profile is given by: 
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where Ec and vc are the composite elastic modulus and Poisson's ratio, 
respectively. 

For a partially bridged crack the distribution of the pressure function 
p(x) is d ivided into two parts: 

unbridged region p(x)= o 00 , for 0 < x <ag ( 2a ) 

bridged region p(x)= o® - c(x), for aq < x <a (3b) 

where an and a are the unbridged crack length and the total crack length, 
respectively. The pressure distribution c(x) is an approximate closure pres- 
sure due to bridging. Equations (1) and (2) can be modified for a partial y 
bridged crack by substituting the appropriate pressures given in equation (3). 

Up to this point both models are identical. The primary difference 
between the two models is the manner by which they account for the closure 
pressure c(x). For the shear lag model, references 3 and 4, c(x) is a func- 
tion of the fiber/matrix interfacial friction shear stress and is proportional 
to the square root of the displacement and given by: 


c(x) 


= 2 


u(x) t vr E c 


R(1 - vf)Em 


1/2 


(4) 


where 

x fiber/matrix interfacial friction stress 
Ef fiber modulus 
Em matrix modulus 
vf fiber volume fraction 

R fiber radius 

By substituting equation (4) into equation (3b), the solution of the displace- 
ment field, equation (2), becomes an integral equation which has to be solved 
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numerically using an iterative scheme with a small damping factor to guarantee 
convergence. After reaching a converged displacement profile, the composite 
SIF for the shear lag model is calculated using equation (1). Finally, the 
effective crack driving force is related to the composite SIF by: 


K Em 

eff " Ec 


(5) 


For the crack closure model, c(x) is equal to the stress carried by the 
fibers in the bridged region averaged over the total bridged area and is given 

b y : / \ 


c(x) 
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( 6 ) 


Equation (6) represents the normal and bending stresses in the bridged fiber 
region and is valid only for a partially bridged crack. The closure load in 
equation (6) is independent of the displacement, hence the SIF for the closure 
model can be solved directly without an iterative procedure. The effective 
crack driving force, K e ff, for the closure model is calculated by substituting 
the results of equation (6) into equations (3) and (1). 


The values of the parameter used for the composite system tested are given 
in table 1. The value of the interface shear stress x depicted in the table 
is a calculated value from experimental indentation tests (ref 3). 


DISCUSSION 

The calculated results of the AK e ff values for both bridging models are given 
in figure 4 as a function of the crack length. Also shown is the AK app for 
the unbridged case. As shown in the figure, the AK e ff values for born 
models decrease with an increase in the crack length. The decrease in the 
AK e f f values corresponds to the observed behavior where the crack growth 
rates decreased and finally crack growth stopped with an increase in the crack 
length. The AK e ff for the shear model is slightly lower than AK e ff for 
the closure model. 

The experimentally obtained composite FCG data was replotted in terms of 
the effective SIF, AK e ff, calculated using both crack bridging models. As 
shown in figure 2, the predicted crack growth rate data for both models are in 
agreement with the data trends exhibited by the unreinforced matrix material. 

An additional check of the accuracy of the bridging models are the numerically 
determined displacements near the crack tip. Both the predicted and observed 
displacements were plotted in figure 5 for a given crack length as a function 
of the distance behind the crack tip. The predicted crack tip displacement 
fields by the two methods are almost identical and are difficult to distinguish 
from each other. Excellent correlation is observed between the measured and 
numerically determined displacement range in the crack tip area. As the dis- 
tance from the crack tip increases some disagreement is observed for both 
mode 1 s . 
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SUMMARY OF RESULTS 


Two fracture mechanics based fiber bridging models were successfully used 
in an attempt to describe experimentally observed fatigue crack growth behavior 
in a SiC/Ti-15-3 composite. The crack closure model and the shear lag model, 
in agreement with the experimentally observed behavior, predicted a decrease 
in the crack growth rates and eventual crack arrest with an increase in the 
crack length. Both models also predicted the near crack tip displacements. 

The advantage of the closure model is that no numerical iteration is required 
to obtain the solution. The success of these models points to the applicabil- 
ity of fracture mechanics to describe the fatigue behavior of certain 
composites . 
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CRACK LENGTH, mm 


TABLE I. - SiC/Ti-15-3 METAL MATRIX PARAMETERS ASSUMED 


o^mflY = 220 MPa 

°°°min - 22 MPa 

R ratio = 0.1 

ag = 1 mm 

w = 5 . 12 mm 

b = 2 .03 mm 

Ef = 427 GPa 

Em = 89 GPa 

Ec = 184 GPa 

R = 72.5 pm 

vf = 0.36 

x = 276 MPa 


CRACK ARREST 


O SHEAR LAG 
■ CLOSURE MODEL 
o PURE MATRIX (TI-15-3) 


O SIDE A 
□ SIDES 













Figure 5. - Crack opening displacements as a function of the distance 
to the crack tip. 
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